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Abstract 

We present a calculation of the fluxes of muons and muon neutrinos from the 
decays of pions and kaons that are themselves the decay products of charmed 
particles produced in the atmosphere by cosmic ray-air collisions. Using the 
perturbative cross section for charm production, these lepton fluxes are two to 
three orders of magnitude smaller than the fluxes from the decays of pions and 
kaons directly produced in cosmic ray-air collisions. Intrinsic charm models 
do not significantly alter our conclusions, nor do models with a charm cross 
section enhanced in the region above an incident cosmic ray energy of 1 TeV. 

I. INTRODUCTION 

The fluxes of atmospheric neutrinos and muons with energies larger than 100 GeV, 
which result from cosmic ray interactions with air molecules in the atmosphere PH, are 
measured in large underground detectors ||. The fluxes are interesting for a variety of 
reasons, including the fact that the leptons are the decay products of particles produced in 
collisions with center of mass energies that may be extremely high. Cosmic rays, comprised 
primarily of protons, have been measured with energies up to E ~ 10 20 eV [|J, many 
orders of magnitude beyond the energies accessible in accelerator laboratories. In addition, 
atmospheric neutrinos and muons are a background for galactic and extra-galactic sources 
of neutrinos || . 

Atmospheric lepton fluxes come from two main sources. The so-called "conventional" 
atmospheric fluxes of muons and neutrinos come from hadronic production of K s, Ki 
and vr^'s in cosmic ray-air collisions, followed by their leptonic, and in the case of kaons, 
semileptonic, decays. At the energies considered here, E > 10 2 GeV, the magnitudes of 
the conventional lepton fluxes are governed by the lifetimes of the mesons. Because of time 
dilation, the decay lengths of these mesons are much larger than the depth of the atmosphere. 
A small fraction decay to leptons, but since they are copiously produced, the conventional 
flux dominates the lepton flux at the lower range of energies considered here. Given an 



incident cosmic ray flux at the top of the atmosphere of 4>cr ~ E a , the flux of leptons 
from K and 7r decays has the energy dependence (j) v ^ ~ E~ a ~ l . Values of a are measured 
with a ~ 2.7-3 [|. 

Charmed particle production and decay make important contributions to the lepton 
fluxes at high energies. Essentially all of the charmed particles decay in the atmosphere 
since r ~ 10~ 12 s, so these contributions are said to contribute to the "prompt" atmospheric 
flux. The energy behavior of the prompt fluxes is one power higher than the conventional 
fluxes for E < 10 7 GeV: <p u ^ ~ E~~ a . Even though charm production is suppressed relative 
to pion and kaon production, the energy behavior of the prompt flux means that the flux 
of leptons from charm decay dominates at high energies. In the recent work of Thunman, 
Ingelman and Gondolo ||, labeled here as TIG, a Monte Carlo model of atmospheric lepton 
production based on PYTHIA |7j was used to calculate the contribution of charmed particle 
semileptonic decays to the atmospheric flux. 

Along with the semileptonic decays of charmed particles, in principle, secondary decays 
such as D — > K — > fiv^ also contribute to the atmospheric lepton fluxes. Essentially all 
charmed meson decays have at least one kaon or pion in the final state, and the charged 
kaons and pions decay mainly to muons and neutrinos. It is the fluxes of these secondary 
muons and neutrinos from charm decay, and the more general topic of the uncertainties in 
the theoretical predictions for the atmospheric lepton fluxes from charm decay, that are the 
topics of this paper. We show that using the perturbative cross section for charm production, 
the secondary lepton fluxes are significant for E < 10 3 GeV relative to the prompt fluxes, 
however they are small compared to the conventional fluxes. Even for strongly enhanced 
charm production, as long as the cross section is consistent with experimental measurements, 
the secondary lepton flux is negligible compared to the conventional flux. 

In Section 2, we outline the standard calculational procedure, an approximate analytic 
method. Details of this procedure are outlined in Refs. 01111. We list our numerical inputs 
in Section 3 and present our results in Section 4. We follow with a summary in Section 5. 



II. CALCULATIONAL PROCEDURE 

The calculations presented here are based on the cascade equations for baryon, meson 
and lepton fluxes @,§,|9] . The cascade equations describe the production of particle j through 
interactions and decays of particle k, and the loss of particle j through absorption and its 
own decay. The energy dependent flux of particles of type j, <pj{E, X), is governed in general 
by the equation 

J2S(k^j) (2.1) 



dX A,- A (dec) 
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where X is the column depth (slant depth) in the atmosphere measured from the top of the 
atmosphere, Aj is the hadronic interaction length and A^- ~ ^ CT . j s the decay length, all 
converted to units of g/cm 2 . The column depth is dependent on the incident angle. For 
vertical fluxes evaluated at sea level, the column depth is X = X = 1300 g/cm 2 . The 
quantity S(k — > j) describes the source of particles of type j from interactions or decays of 
particles of type k. For neutrinos and muons at energies above 100 GeV, since the interaction 
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and decay lengths are essentially infinite, only the source terms appear on the right hand 
side of equation (2.1). The source terms have the form: 



oo 



(,2) 

where A^(^) is an interaction length or decay length, depending on the source. The quantity 
dnk^j/dE describes the energy distribution of the produced particle j, and (Nj) is the 
particle multiplicity. For decays, the value of (Nj) is unity, except for pions. For the 
production of charm, (N c ) = 2. 

In principle, the cascade equations are a coupled set of equations. We make several 
simplifying assumptions so that we can solve the equations analytically. First, we decouple 
the meson equations, with the exception of charm decay to pions and kaons. In principle, 
there are interaction terms such as S(n — ► K), however, these are small compared to S(p — > 
K), so we neglect the former. With this assumption, we can solve for conventional 7r and 
K fluxes separately, then evaluate their contributions separately to the neutrino and muon 
fluxes. The same holds true for charm production. We follow the produced charm in their 
decays to leptons, pions and kaons. 

To further simplify the calculation, we represent the cosmic ray flux by the proton flux. 
The cosmic ray composition is dominantly protons [10[], so the treatment of cosmic rays as 



protons is reasonable. This approximation was used by TIG in Ref. 0. A feature that 
allows us to approximately solve the cascade equations is to use a factorization of the energy 
dependence in the fluxes. By writing 

<f>j(E,X)=fj(E,X)E-^ , (2.3) 

and assuming that fj(E,X) is weakly dependent on energy, approximate solutions to the 
cascade equations are straightforward and lead to the power law energy behavior described 
in the introduction. Thunman, Ingelman and Gondolo have demonstrated this is a good 
approximation by comparing the approximate solutions to Monte Carlo results ||. This is 
especially useful because the source terms factorize to 
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S(k - j) = Jk \ ; m Z kj ((l k , E) (2.4) 



where Z k j((5 k ,E) depends only weakly on energy. The calculation of the Z-moments is the 
essential ingredient in addition to A to describe the lepton fluxes. The interaction moments 
include information about multiplicity, the magnitude and energy dependence of the relevant 
cross section and the energy distribution of the emerging particle. To simplify notation, we 
omit the energy argument in Z's, although we do include energy dependent Z's in our 
evaluation of the fluxes. 

With the assumptions above, the solutions for the hadron fluxes are straightforward. For 
cosmic ray protons, \ p dec?> is infinite, and the solution for the flux is 

<j> p (X, E) = f p (X)E- a = f p (0)e- x ^?E- a (2.5) 

for A*; = Afc/(1 — Zkk). For hadrons other than the proton, the solutions to the factorized 
cascade equations have two different forms, depending on whether one is in the high or 
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low energy regime. The energy regime is determined by whether or not the decay length 
of the particle is large compared to the column depth. The cascade equations are solved 
approximately in these two regimes, then interpolation between high and low is done with 



Idowjhigh 



, 3 yj .. . . (2.6) 
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Explicit solutions appear, for example, in Ref. |J. The fluxes from each source are evaluated 
this way, then summed to get the total flux. 

The lepton fluxes come from hadron decays, where the hadron fluxes are approximated 
via Eq. (2.6). In the high and low energy regimes, the approximate lepton fluxes can 
be written in simple form. For definiteness, we write the flux for neutrinos. The flux for 
muons is in exactly the same form. At high energies, the (direct) flux from source terms 
p — » j ^leptons is 

^ h (X,E) = fi ,hi9h (E, X)E~ a ~ 1 = z vM z M)H^IK) fM E - a (27) 
v K ' v y J 1-Z pp (a) 1 - VA,- Af c) p(X) y ' 

where p(X) is the density of air at column depth X. The energy dependence of A^ 6 ^ makes 
the flux of leptons from high energy hadron decays suppressed by one power of energy 
relative to the proton flux. For low energies, the hadron flux itself is proportional to X^f^ 
and the decay length cancels. The energy dependence of the lepton flux tracks that of the 
primary cosmic ray proton flux. From the direct decays of low energy hadrons j , 

<% 1 ™(X, E) = f^(E, X)E~ a = y|f^ / F (0)^ , (2.8) 

1 - Z pp [a) 

As indicated earlier, pions and kaons have relatively long lifetimes, so for E > 10 2 GeV, 
they are "high energy" hadrons resulting in a lepton flux proportional to E~ a ~ x . Charmed 
hadrons, for most of the energy range E ~ 10 2 — 10 8 GeV, are low energy particles, so the 
resulting lepton flux has the low energy form proportional to E~ a . 

The extensions of these equations for secondary leptons from p — >charm^ j — deptons 
is straightforward. The substitution 

Z pj (a) -> Z pkc (a)Z kcj (p kc ) (2.9) 

is required for hadrons j that come from charmed particle (k c ) decays. The calculation of the 
moments for direct processes have be done with a variety of assumptions ||[2[ . It remains to 
indicate the inputs for interaction lengths and Z-moments and to describe our calculation 

of Z kcj (f3 kr ). 



III. INPUTS 



Eqs. (2.7) and (2.8) govern the lepton fluxes, with the additional substitution of Eq. 
(2.9) for the secondary fluxes. We use TIG parameters where they are available 0. TIG 
have shown that for conventional charm production, D ± , D° and D° decays dominate the 
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sources of prompt leptons for most of the energies between 10 2 — 10 8 GeV. Consequently, we 
consider charm contributions only via D ± , D° and D° decays. Henceforth, we will only refer 
to particles, however, all of our calculations include both particles and antiparticles. The 
inputs fall into three categories: the incident proton flux (we assume that cosmic rays are 
protons), interaction Z-moments and interaction lengths, and decay Z-moments and decay 
lengths. 

The incident flux of protons at the top of the atmosphere, following TIG, is taken to be 
isotropic, with 

p (E,X = O)[cm-V 1 sr- 1 GeV- 1 ] = 1.7 (E/GeV)~ 2 - 7 E < E (3.1) 

174 (E/GeV)- 3 E > E , 

with E = 5 • 10 6 GeV. 

The interaction lengths for the protons, pions and kaons require the total scattering cross 
sections with isoscalar nucleons. We use the Donnachie-Landshoff parameterization [[11|] 

a tot = Xs* + Ys-v (3.2) 



and the values of the parameters from the Particle Data Book [12] for incident protons, n + 
and K + . We use the K + cross sections for Kl, D° and D + as well. The cross sections are 
then rescaled with a constant factor of A 2 / 3 , where A = 14.5 is the average atomic number 
of air nuclear targets. Interaction Z-moments have been presented as a function of energy 
in Ref. || for nucleons, pions, kaons and D mesons. We use these Z-moments to evaluate 
Eqs. (2.7) and (2.8). 

For particle lifetimes and branching fractions, we use the Particle Data Book values. 
The two-body decay moments for pion and kaon decays are straightforward to calculate. 
Formulae appear in Ref. M, and numerical values, depending on the energy behavior of 
the decaying mesons, appear in Refs. || and @. For semileptonic and non-leptonic decays, 
the calculations are less obvious. These must be done for the decay moments for charm 
to pions and kaons, muons and neutrinos. An accurate evaluation of Zk c j(/3 C ) is difficult 
because the hadronic decay distributions cannot be calculated analytically. Data exist on 
distributions in charmed particle decays, for example, the K ± and K° momentum spectra 



for D°, D + and Df, where the .D's are all produced approximately at rest |T3[. The data are 



not immediately translatable into decay Z-moments because one must boost to the frame 
where the kinetic energies of the .D's are large. 

For the simplest case of semileptonic charm decays, TIG distribute lepton momenta 
according to a weak matrix element 

I M | 2= (Pd ■ P^){Pu ■ Ph) (3.3) 

where h is the final state hadron system. They find that the decay moments to muons and 
neutrinos are equal to within 15% for D mesons, and that the moments decrease by a factor 
of ~ 3 as the .D's make the transition from the low energy to high energy regime. 

The matrix element in Eq. (3.3) does not account for form factor modifications of 
the hadron weak vertex. To look at the sensitivity of the decay moment to the energy 
distribution, we have calculated the decay Z D ± U using three-body phase space to dictate the 
energy distribution of the neutrino. Our results are shown in Fig. 1 by the points. The 



5 



energy dependence comes from the changing energy behavior of the flux. The upper 
dashed line is the TIG result with (3d = 1.7, appropriate for low energy .D's, for example, 
at E = 10 2 GeV. The lower dashed line is with (3d = 3, appropriate for high energy .D's, 
namely at 10 8 GeV. Over the range of energies between 10 2 — 10 8 GeV, there is a factor 
of ~ 2 — 2.4 between the phase space result and the relevant TIG numbers. This factor of 
2 — 2.4 gives a qualitative estimate of the uncertainty in the decay moments. We use the 
TIG semileptonic decay moments, interpolating between 10 2 and 10 8 GeV as a function of 
\n(E/GeV) using a constant (negative) slope. 
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Fig. 1. The points indicate the decay moment Z D ± V evaluated using three-body phase space 
to distribute the neutrino energy. The upper dashed line is the TIG result using (3 D = 1.7. 
The lower dashed line has (3d — 3 0. 



Our procedure for the hadronic decay moments of the .D's is to use the same energy 
dependence and magnitude as the semileptonic D moments as a function of energy, with a 
correction to account for the different branching fractions and particle multiplicities of the 
hadrons. For example, 

(N K ±) BjD* ^ K± + anything) 
Zd±k± = Zd±v BW - v) ( } 

The Zk c j(f3 c = 1.7)-moments and the branching fractions and multiplicities used are 
shown in Table I. In the conversion to the appropriate Z-moments, the particle data book 
values for the hadronic branching fractions [|12|] were used. In addition, we use Mark III 



results on the topological branching fractions The branching fractions for .D's into 

kaons is in the Particle Data Book. For pions, we take the branching fraction into pions 
equal to the fraction of all decays which are not semileptonic. The kaon multiplicities for 
D° and D + decays are assumed to be (N K ) = 1. To estimate (N w ±), we use the charged 
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particle topological branching fractions. Mark III data Jl3| indicate that the average charged 
particle multiplicities for D° and D + range between 2.4 and 2.6. We estimate that the 
average hadronic charged particle multiplicities are between 2.1-2 A. Assuming that either 
K , K° or K° accompanies all charmed meson decays, and that semileptonic decays have 
no pions, we estimate that the average number of charged pions in the decays are between 
1.5-1.8, as indicated in Table I. 



IV. RESULTS 

The parameters of Table I complete what we need to evaluate the lepton fluxes. In Fig- 
ures 2 and 3, we show our resulting muon and muon neutrino (particle plus anti-particle) 
fluxes scaled by a factor of E 3 , evaluated in the vertical direction. Our conventional, prompt 
and secondary fluxes are shown by the solid lines. The dashed lines are the TIG parameter- 
izations which fit their Monte Carlo results for the conventional and prompt fluxes M. 
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Fig. 2. The solid lines show the vertical muon plus antimuon flux, scaled by E 3 , as a function 
on muon energy E for conventional, prompt and secondary decay sources. The dashed lines 
show the TIG parameterization of the conventional and prompt fluxes. 

Before discussing the secondary fluxes, we comment on the discrepancies between the TIG 
parameterization and the prompt lepton fluxes above E = 10 5 GeV. Part of the discrepancy 
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comes from our ignoring D s and A c contributions. TIG have shown that the A c , in particular, 
makes significant contributions at the higher energies. In addition, while we have attempted 
to use the same parameters as TIG, not all of them are specified in their paper. Finally, 
their fit was made to Monte Carlo results, which are not exactly reproduced by this semi- 
analytic method using Z-moments. In any case, our focus is on secondary leptons in an 
energy regime where our results are in reasonable agreement with TIG's parameterizations. 

We note that the secondary muon flux is larger than the prompt muon flux for E < 10 3 
GeV. The secondary muon neutrino flux exceeds the prompt neutrino flux only for E < 200 
GeV. The flux of secondary muons (and conventional muons) is larger than the corresponding 
neutrino flux because muons in pion decay carry the bulk of the pion's energy. For fixed 
energy of the lepton, the parent pion energy is lower for muons than for muon neutrinos. 
This is a significant feature when the parent pion flux is decreasing rapidly with energy. 




Fig. 3. The vertical muon neutrino plus antineutrino fluxes, as in Fig. 2. 



The secondary p — > charm — > n, K ^lepton chain is suppressed by a factor of 

Rj = EfccZpfcAj (4.1) 
Z P j 

relative to the conventional p — > ir,K — depton chain. For TIG's calculation of charm 
production based on PYTHIA, Z p k c ~ 10~ 3 Z p7T ~ 10~ 2 Z pK . Our secondary fluxes are 
suppressed by a factor of ~ 1/(2.4 x 10 3 ) for muons and ~ 1/(1-1 x 10 3 ) for muon neutrinos 
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relative to the conventional fluxes. The errors due to our approximations in the calculation 
of the conventional fluxes are more significant than the secondary flux as long as the charm 
interaction and decay moments are reliably known. We turn now to the question of whether 
or not the factor Rj can be increased such that the secondary flux is a significant contribution 
to the total flux below lepton energies of 1 TeV. 

Increasing Rj may be achieved by increasing the hadronic decay moments and the inter- 
action moments associated with the production of .D's. We will proceed by estimating an 
overall enhancement factor K with 



K = K d K a K x , (4.2) 

the product of enhancement in decay (d), cross section (a) and in the charm energy dis- 
tribution (x = E c /E p ), which may multiply Rj. We'll assume that these are universal to 
all Rj. For the limited energy range where secondary lepton fluxes are larger than prompt 
fluxes: E = 10 2 — 10 3 GeV, we take K independent of energy. A full numerical study of the 
uncertainty in the charm production moments as a function of energy is in progress WA . 

As described in Sec. Ill, we have evaluated the decay moments by rescaling the V — A 
Zj) v moments according to multiplicities and branching fractions. By using three body 
phase space, the decay moments are enhanced by a factor of 2.4 at E = 10 2 GeV. We take 
K d = 2.4. 

The cross section at next-to-leading order (NLO) in QCD is not very well specified 
because of the low mass of the charm quark. A variation of the charm quark mass between 
1.2 and 1.8 GeV, and the renormalization scale between m c /2 and 2m c in the NLO cross 



section gives a range of cross sections that vary by more than two orders of magnitude |15 
By comparing the TIG cross section with fixed target cross section results, a factor K a = 2 
is not unreasonable. 

The x distribution is an important input into the calculation of the interaction Z- 
moments because of the steeply falling fluxes. The quantity x = E c f E p defined in the 
fixed target frame is approximately equal to | xp |= 2 | pjj m | /y/s, in terms of the longi- 
tudinal momentum of the charmed quark or antiquark in the hadron center-of-momentum 
frame. The theoretical (perturbative QCD) distributions in xj? are typically softer than 
the measured distributions, which would tend to lower the predicted Z-moment. The xf 
distributions for xf > are parameterized by da/dxF ~ (n + 1)(1 — x^) n . We make the 
scaling assumption so that the interaction Z-moments are proportional to 

1 dxx 17 A n (l -x) n (4.3) 

for the proton flux proportional to E~ 2 - 7 and normalization factor A n = n + 1. For low 
mass charmed quarks (m c = 1.2 GeV), the NLO calculations of the xp distributions, fit to 
a (1 — xf) h form, give n ~ 6 — 9.5 for E = 10 2 — 10 3 GeV [|15|. Fixed target data yield 



n ~ 4.9 — 8.6 | 16fl , some with large error bars. To estimate K X} we take the NA32 value 



for n, n — 5.5 [T7| for a typical experimental value, and n = 7.5 as representative of the 



perturbative value, so 



Jo 1 dxx L7 6.5(l 



x 



,5.5 



K x = ^ ^ '— = 1.5 . (4.4) 

Jo 1 dxx L7 8.5(l-x) 7 - 5 V 1 



9 



Taken together, the overall enhancement factor based on modifications to the TIG charm 
production parameters is on the order of K ~ 7. 

Intrinsic charm, in principle, may boost the charm production rate. Intrinsic charm 
models have a charm (and anti-charm) component in the proton at a low scale Qq, before 



QCD evolution has generated cc pairs by gluon splitting pq| . Estimates of the Z-moments 
for several models of intrinsic charm were made by TIG |J , and maximum Z-moments of a 
few times 10 -3 were found. These are at just the level of the perturbative charm Z-moments, 
so again, may account for at most a factor of two. 

Finally, we turn to a model which does not rely exclusively on the perturbative parame- 
ters presented by TIG and energy independent modifications. We use a charm cross section 
at high energies suggested by Zas, Halzen and Vazquez Jl9| where 



a c5 = 0.1a tot , (4.5) 

where Otot appears in Eq. (3.2). The form of the cross section in Eq. (4.7) seriously 
overestimates the cross section in the measured regime, E = 10 2 — 10 3 GeV. We take the 
hadronic production cross section to be 

where er^ is the next-to- leading order perturbative cross section with m c = 1.3 GeV, 
evaluated at the renormalization and factorization scales set to \x = m c , for E < 10 3 GeV, 



using the CTEQ3 parton distribution functions p0[ . This is a slight overestimate of the 
data below E = 10 3 GeV. Above 10 3 GeV, we take a power law 



a™ = 1.3 • 10 5 (£/GeV) a8b5 pb , (4.7) 

which is the power law extrapolation of the perturbative cross section at E = 10 3 GeV. 

By E ~ 10 6 GeV, <r c5 ~ 0.1a tot . This assumption yields significantly higher prompt 
fluxes at high energies. Indeed, Gonzalez-Garcia, Halzen, Vazques and Zas |21| have shown 



that such large charm cross sections above energies of 10 GeV are problematic when one 



looks at Akeno horizontal air shower data [22]. For our purposes here, we are interested in 



the cross section only as it affects the interaction Z-moments below 1-10 TeV. For the x 
distribution, we assume a scaling form, da/dx ~ 5(1 — x) 4 , which is is the lower range of 
consistency with measured values p3| . Using these inputs, our results for the prompt and 
secondary muon fluxes from charm decay are shown with the dashed lines in Fig. 4. For 
comparison, the solid lines indicate the TIG results and our secondary flux calculation based 
on TIG parameters, as in Fig. 2. The difference between the fluxes at E = 10 2 GeV are 
due to differences in total cross sections and x-distributions relative to PYTHIA generated 
distributions. 
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Fig. 4. The vertical muon prompt and secondary muon fluxes, scaled by E 3 , using the 
conventional TIG parameters (solid line) and using Eq. (4.8) and da/dx ~ 5(1 — x) 4 
(dashed line). 



It is clear from Fig. 4 that at energies above a few TeV, the high energy behavior of the 
cross section and x-distribution have important implications for the prompt neutrino flux. 
Recall that at these energies, the prompt muon flux equals the prompt muon neutrino flux, 
and also equals the prompt electron neutrino flux. The prompt fluxes are isotropic. For the 
enhanced prompt flux of Fig. 4, the crossover from the vertical conventional-dominated to 
prompt-dominated muon neutrino fluxes occurs around E ~ 3 • 10 4 GeV. This is in contrast 
to TIG's result of a crossover close to E ~ 10 6 GeV. However, there is no significant effect 
for E < 10 3 GeV using the enhanced cross sections, given our overall uncertainty of a factor 
of K ~ 10. Indeed, below E = 10 4 GeV, a factor of 10 uncertainty in the lepton fluxes from 
charm decay is a reasonable estimate. 

Finally, we remark that while we have focussed on secondary leptons in the energy 
interval 100-1000 GeV, it is interesting to note that the dashed line in Fig. 4, if valid at 
E = 10 4 GeV, has important implications for the electron neutrino flux. At E — 10 4 GeV, 
the prompt muon neutrino flux is still an order of magnitude below the vertical conventional 
flux, however, the prompt electron neutrino flux is larger than the vertical conventional flux. 
This can be seen by considering the ratio R = <pu^/4>u e as a function of zenith angle 9 and 
energy. Using the Bartol group's calculation of the conventional atmospheric lepton fluxes 
24J, neglecting any charm contribution, the ratios are: 



Rno C (E„ = 10 4 GeV, cos 9 = 1) 



32 



(4.8) 
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R no C (E U = 10 4 GeV, cos# = 0) = 29 



However, when one adds in the prompt flux shown by the dashed line in Fig. 4, the same 
ratios yield 

R mthc {E u = ltf GeV, cos = 1) =7 (4.9) 
Rwm c{E u = 10 4 GeV, costf = 0) = 18 . 



Measuring R is a difficult task, especially at 10 TeV, because of the difficulty of measuring 
the electron neutrino flux. Air shower arrays are, in principle, sensitive to electron neutrino 
induced events, however, acceptances are low compared with higher energies |25j . In under- 



ground experiments, one does not have the advantage of the long muon range that increases 
the effective volume of the detector when the process is v e N charged-current interactions. 



V. CONCLUSIONS 

We have focused here on corrections to the atmospheric lepton fluxes from charm decays 
in the energy range E = 10 2 — 10 3 GeV, where secondary neutrinos from the D — ► K,tt — > 
decay chain contribute. The secondary flux corrections are significant for the muon 
flux compared to the prompt muon flux below ~ 10 3 GeV, and for the muon neutrino fluxes 
below ~ 200 GeV. Since the branching fraction for K and 7r ± decays to electron neutrinos 
is small, secondary corrections to the atmospheric electron neutrino fluxes is negligible. 

The prompt and secondary lepton fluxes are tied to the charm cross section and energy 
distribution of the emerging charmed particles. We have demonstrated here that the exper- 
imental constraints on the charm production cross section and energy distribution limit the 
extent to which leptons from the decays of charm can contribute to the overall atmospheric 
lepton fluxes. Taking the parameters of TIG in Ref. [2] as a starting point, we have shown 
that the uncertainties in the calculation of the lepton fluxes from charm decay may result 
in an enhancement by as much as a factor of 10 in the energy range of 10 2 — 10 3 GeV. 
An additional two orders of magnitude are required to make the secondary lepton fluxes 
comparable to the conventional ones at E = 100 GeV. At higher energies, where one is less 
constrained by experiments, the high energy behavior of the charm cross section and energy 
distribution significantly influences the prompt flux, a topic of further investigation |TJ . 
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TABLES 



TABLE I. Charmed particle decay moments. 
Decaying Particle (fc c ) i B (N^ Z kci {(3 c = 1.7) 



D+ 




0.30 


1 


0.032 




K L 


0.30 


1 


0.032 




7T ± 


0.66 


1.5 


0.10 






0.17 


1 


0.018 


D° 


K± 


0.56 


1 


0.069 




K L 


0.21 


1 


0.026 






0.85 


1.8 


0.19 






0.068 


1 


0.0084 
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